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Discussing the general case of a hard partonic production process, we show that the notion of 
parton energy loss is not always sufflcient to fully address medium-induced gluon radiation. The 
broader notion of gluon radiation associated to a hard process has to be used, in particular when 
initial and final state radiation amplitudes interfere, making the medium-induced radiated energy 
different from the energy loss of any well-identified parton. Our arguments are first presented 
in an abelian QED model, and then applied to large-a;^^ quarkonium hadroproduction. In this 
CEise, we show that the medium-induced radiated energy is qualitatively similar (but not identical) 
to the radiative energy loss of an "asymptotic massive parton" undergoing transverse momentum 
broadening when travelling through the nucleus. In particular, it scales as the incoming parton 
energy, which suggests to reconsider gluon radiation as a possible explanation of large quarkonium 
suppression in p-A collisions. We expect a similar effect in open heavy-flavour and possibly light- 
hadron hadroproduction at large , depending on the precise definition of the nuclear suppression 
factor in the latter case. 



PACS numbers: 12.38.Bx, 13.85.-t, 25.75.-q 



INTRODUCTION AND MAIN RESULTS 



High-energy collisions involving atomic nuclei give a 
unique opportunity to study parton propagation in nu- 
clear matter. Deep inelastic scattering (DIS) on nuclei 
and the Drell-Yan (DY) process in proton-nucleus colli- 
sions are usually considered to be good probes of parton 
propagation in cold (confined) nuclear matter, whereas 
large-p^ hadron or jet production in heavy-ion collisions 
should be sensitive to parton propagation in hot quark- 
gluon plasma (QGP). It seems clear that in order to 
provide a sensible interpretation of the phenomenon of 
jet-quenching, observed at RHIC jT| and recently at the 
LHC [2j, and considered as a prominent QGP signal, a 
good theoretical understanding of parton propagation in 
cold nuclear matter is required. 

An energetic parton travelling in a large nuclear 
medium undergoes multiple elastic scatterings, which in- 
duce gluon radiation. Usually, the amount of such gluon 
radiation is referred to as the (radiative) energy loss of 
the fast parton, and the notion of parton energy loss has 
been widely used in phenomenological studies of nuclear 
effects. For instance, the suppression of production rates 
in A-A or p-A compared to p-p collisions (after an ad- 
equate normalization) observed for various processes in 
some kinematical regions (DIS at large 2, large- a;^ DY 
production in p-A, large-p^ hadron production in A-A) 
has been attributed, at least partly, to parton energy 
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loss [3]. Intuitively, due to parton energy loss, a hard 
QCD process probes the incoming parton distribution 
functions (PDF) at higher x where they are suppressed, 
leading to nuclear suppression. 

However, there exist some hard processes where the 
medium-induced associated radiation cannot be strictly 
identified with the energy loss of a well-defined parton. 
Such processes are those where a color charge (which may 
be a composite object) is produced nearly collinearly with 
one of the incoming partons (in the rest frame of the 
medium, see section |II A), making the gluon emission 



amplitudes before and after the hard production vertex 
coherent. Quite intuitively, the presence of interference 
precludes the identification between associated radiation 
and radiation off a well-defined parton. The notion of 
radiated energy associated to a hard process is thus more 
general than the notion of parton energy loss. The aim 
of our study is to illustrate the latter statements, using 
the specific example of large- a;^ quarkonium hadropro- 
duction in high-energy p-p and p-A collisions. 

Our main result can be summarized as follows. The 
medium-induced gluon radiation associated to large- a:;^ 
quarkonium hadroproduction arises from large gluon for- 
mation times, scales as the incoming parton energy 
(which at large becomes commensurate with the 
quarkonium energy), and cannot be identified with the 
energy loss of a well-defined parton. However, it is qual- 
itatively similar to the Bethe-Heitler energy loss of an 
"asymptotic" massive parton. A similar effect should 
arise in all processes where a fast color charge is pro- 
duced nearly collinearly with a fast incoming parton. 
For instance, the effect is expected in large-x^ open 
heavy flavour production, but not in Drell-Yan produc- 
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tion, where the final energetic particle carries no color 
charge (see section IV for a more detailed discussion). 



frame) is large compared to its transverse mass 



In order to motivate our study, let us first shortly re- 
view the notion of parton energy loss. The latter has been 
extensively studied (see Ref. |4] for a heuristic review). 
It is of course of crucial importance to know the correct 
parametric dependence of the mean parton energy loss 
Ai5, which in general depends on the parton properties 
(in particular its energy E and mass M), those of the 
medium (its size L and density and/or temperature), and 
on the coupling as . But as emphasized in Ref. [Ij , the ra- 
diative loss of an energetic charged particle also depends 
on the way the particle is produced. When discussing 
radiative energy loss, two physical situations have to be 
distinguished: 

(i) the energetic charge is suddenly accelerated in the 
medium, at a given initial time t = (as in DIS 
for instance). Since a suddenly accelerated particle 
radiates even in vacuum, in this case the relevant 
quantity to address nuclear suppression is not the 
total energy loss AE suffered by the parton, but the 
additional, medium-induced loss Ai5i„d — AE'^cd — 
AE'vac occurring in the medium when compared to 
vacuum; 

(ii) the energetic charge is an asymptotic particle pre- 
pared at t = — oo and detected at t — -foo, travel- 
ling in the medium between t = and t = L. In 
this situation, the medium-induced loss coincides 
with the total loss. 
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the QQ pair is ultra-relativistic and produced nearly 
coUinearly to the projectile parton. Moreover, at large 
enough p^^ , the QQ pair remains compact (of size 
1/M±) and in a color octet state for a time ioctet ^ L |42| . 
as a result of time dilation in the nucleus rest frame. Un- 
der those conditions, the hard production process looks 
like gluon scattering at small angle in the target rest 
frame, except that the outgoing "gluon" - the octet 
QQ pair - is massive. This intuitively explains why the 
medium-induced radiated energy associated to large- a;^ 
quarkonium hadroproduction is qualitatively similar to 
the medium-induced energy loss of an "asymptotic mas- 
sive parton" crossing the target nucleus (see section III I . 



In particular, its scaling as the quarkonium energy might 
explain the strong nuclear suppression of quarkonium 
hadroproduction observed at large Xp . 

The suppression of J/^ production in hadron-nucleus 
collisions at large has been observed by several ex- 
periments, at different collision energies y/s~^ ~ 20- 
40 GeV/nucleon and in various nuclear targets [IHS]. 
Remarkably, its magnitude is similar to that of open 
charm measured by E866/NuSea through single muon 
production |9j , as well as that of light hadrons measured 
by NA49 [lO] at SPS and BRAHMS [H] at RHIC; in 
contrast, it proved much stronger than observed in the 
DY channel [H] as shown in Fig. [l] where J/ip and DY 
E866/NuSea data are plotted as a function of Xp. 



The situation (ii) is more natural in QED, where 
asymptotic charges do exist. In QCD, the only way 
to "see" a color charge is to resolve it inside a color- 
less hadron, via some hard partonic subprocess. Even 
in the absence of a nuclear medium, say in e-p or p- 
p collisions, the hard process inevitably produces gluon 
radiation. Thus, the "medium-induced prescription" is 
required, and at first sight, it is the situation (i) which is 
generic in QCD. For example, in DIS (DY) production, 
the medium-induced associated radiation can be iden- 
tified with the medium-induced energy loss of a quark 
created (annihilated) at t = 0, obtained by subtracting 
the loss associated to the hard process in e-p (p-p) from 
that in e-A (p-A) scattering. Large-p^ hadron produc- 
tion around mid-rapidity in A~A collisions is also sensi- 
tive to the medium-induced loss of a parton created at 
t = 0. (See section II B for a discussion of this case.) 



However, as mentioned in the beginning, for some hard 
processes the associated medium-induced radiation does 
not correspond to the energy loss of a well-defined parton 
- in either situation (i) or (ii). To illustrate this possi- 
bility, let us consider large- a;^ quarkonium hadroproduc- 
tion, and assume for simplicity that the heavy quark- 
antiquark QQ pair is dominantly produced through the 
gluon fusion channel, gg QQ [H]. When the longitu- 
dinal momentum p.. of the QQ pair (in the nucleus rest 
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FIG. 1: Atomic mass number dependence {a{p-A) oc A") of 
J/^p (circles) and Drell-Yan (squares) production (4 < M < 
8.4 GeV) measured in p-A collisions at y/s^ = 38.8 GeV as 
a function of x^. E866/NuSea data from [SI 112). 

Various explanations have been proposed to account 
for these measurements. The comparison of NA3 [S], 
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E772/E866 [ilH] as well as PHENIX [H] measurements 
allowed to rule out a scaling of J/t/i suppression in the 
target parton momentum fraction X2 , indicating a viola- 
tion of factorization in charm hadroproduction [14] . This 
observation consequently ruled out the nuclear modifica- 
tions of parton densities in the target as being a dom- 
inant effect. On the contrary, the observed scaling in 
either or xi supported the qualitative predictions of 
the intrinsic charm picture in which the soft scattering of 
higher-Fock states (e.g. \uudcc)) on the nucleus leads to 
A^/-^ scaling {a — 2/3) at large x^ [13]. Although intrin- 
sic charm might contribute significantly to the observed 
J/ifj suppression even at rather small Xp, its sole effect 
cannot explain the data given the present constraints on 
the amount of charm in the proton |16| . The inelastic 
interaction of the cc pair with nuclear matter seems also 
disfavoured, since it would require unrealistically large 
cross sections ;17]; also, J/ip nuclear absorption models 
would naturally lead to X2 scaling [181 unlike what is ob- 
served experimentally. 

It has also been suggested by Gavin and Milana (GM) 
that parton energy loss in nuclear matter might be the 
dominant effect responsible for both DY and J/ip sup- 
pression [15]. We believe that the present paper should 
resolve some confusion concerning their explanation. In 
the GM model, the energy loss is assumed to scale with 
the incoming parton energy (without real justification), 
AE (X E, leading to a shift Axi cx Xi in the projectile 
PDF. This is a key assumption in order to reproduce the 
qualitative trend of the data. If, on the contrary, energy 
loss was independent of the parton energy, the shift in 
the momentum fraction would not depend on xi, thus 
leading to a less steep J/t/j suppression as a function of 
xi, unlike the trend of the data [12]. The stronger J/ip 
suppression, compared to that of DY, arises from several 
effects in GM [19j. Since J/tp production proceeds mostly 
via gluon fusion (at not too large Xp), the typical energy 
loss is expected to be larger than for quarks by a factor 
Nc/Cf = 9/4. Another effect comes from the large-x 
PDF being steeper for gluons than for quarks. Finally, 
the produced cc pair produced at large Xp travels as a 
color octet through the nucleus, thus losing energy as a 
gluon provided the multiple collisions in the final state do 
not resolve the cc pair. The latter effect appears to be 
quantitatively crucial as it accounts for more than half 
of the J/ijj suppression. In summary, in the GM model, 
J/^jj suppression is due to the initial (quark or gluon) and 
final (color octet cc pair) in-medium energy losses added 
incoherently. 

With this in mind, it seems clear that the "parton en- 
ergy loss" implied in the GM model refers to the medium- 
induced energy loss of a color charge created (or annihi- 
lated) at the hard production time t — 0, i.e., to the sit- 
uation (i) defined previously. As we will now recall, such 
parton energy loss does not scale with the parton energy 
E when E oo (all other scales being fixed), which ap- 
parently rules out the Gavin-Milana model. Consider a 
parton created at some given time t = 0. As already men- 



tioned, the initial acceleration produces radiation since 
the initially bare parton progressively builds its proper 
field, a process which goes along with associated radi- 
ation. This happens even when the parton is created 
in vacuum. In the presence of a nuclear medium, the 
medium-induced radiation arises due to the rescatterings 
off target partons shaking off the components of the par- 
ton proper field just built. Over the distance L, only the 
components of the proper field which are already formed 
can be shaken off, hence the medium-induced radiation 
must satisfy the constraint 
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<L. 



(1.2) 



where uj and k± are the radiated gluon energy and trans- 
verse momentum, ~ fcj^/a; <C 1 the emission angle, 
and tf the gluon formation time. The constraint (1.2) 
yields a bound on parton energy loss, as first derived by 
Brodsky and Hoyer 20 , AE ^ u < L{k\). As noted 
in Ref. [20., the latter bound is a direct consequence of 
the uncertainty principle, ALApz > 1, implying that 
some radiation can be released over the length L provided 
the emission process involves a large enough longitudinal 
momentum transfer, > A p^ > 1/AL > 1/L. With 
Pz ~ kj_/{2u}), one recovers (1.2). We stress that the 
Brodsky-Hoyer bound applies to situation (i), i.e., to the 
medium-induced energy loss of a parton created (or an- 
nihilated) in the medium. In this situation, the medium- 
induced radiation probes the medium size L, as stated by 
(1.2). In other words, large formation times tf ^ L 
cancel out in the medium-induced loss. Those state- 
ments were later confirmed by explicit calculations of the 
medium-induced radiative loss of a parton created in a 
medium, showing that AE cx L^E'^ when E ^ oo [21]. 

As legitimately claimed by Brodsky and Hoyer [20) . 
the bound on parton energy loss seems to rule out the 
Gavin-Milana "energy loss explanation" of J/i/j nuclear 
suppression (which uses ad hoc parton energy losses cor- 
responding to situation (i) and nevertheless scaling as E) . 
However, we believe this conclusion relies on a misinter- 
pretation of the physics at work. As we already empha- 
sized, large-xj? quarkonium production is a process where 
the emission amplitudes off the nearly coUinear incoming 
and outgoing color charges strongly interfere. Thus, the 
gluon radiation spectrum associated to the hard process 
cannot be identified with the radiation spectrum off a 
well-defined parton (or with the incoherent sum of such 
spectra), contrary to what is assumed in Ref. 1 1£| . 

This is best illustrated by our calculation of sec- 
tion II C where large-x^^ quarkonium production is mod- 



elled by a simple QED process, where an energetic muon 
of mass M is produced in the hard scattering of an in- 
coming electron of mass m, see Fig. [7j The (photon) 
radiation spectrum associated to the hard process in a 
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p-p collision is given by (see Fig. [s] and (2.25)) 
dl 
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(1.3) 



where the last term stands for the interference alluded 
to above, 6's|pp = q±/E denotes the angle of the fi- 
nal muon with respect to the incoming electron, and 
6m = M/E, = m/E. The radiation spectrum is simi- 
lar to the Bethe-Heitler spectrum of an asymptotic charge 
of mass M undergoing a scattering of angle 9s — 0s\pp, 



see (2.361, the latter being approximately obtained by 
setting m = M in (|1.3[). Although the associated spec- 



trum (1.3) and Bethe-Heitler spectrum are not strictly 
identical, they share the same property which is essen- 
tial to our discussion. They both arise from large photon 
formation times, tf ^ l/{u!0'^) ^ L, resulting in an inte- 
grated spectrum scaling as E. 

Obviously, the same is true for the spectrum associated 



to the hard process in a p-A collision, obtained from ( 1.3 ) 
by the substitution 6*5! 



0,|pA (see (2.26)) 
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(1.4) 



Our main observation is that the medium- induced spec- 



trum, obtained by subtracting (1.3 ) from ( 1.4 1, also arises 
from large formation times tf ^ L, and moreover solely 
from the interference terms. 
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(1.5) 



where A^^ — ^s|pA ^ ^^Ipp' a consequence, the 
medium-induced radiated energy associated to the large- 
production process scales as the energy E of the 
incoming charge, similarly to Bethe-Heitler radiation. 
This conclusion trivially generalizes to the QCD case of 
quarkonium hadroproduction (see section HI). Thus, the 
Gavin-Milana assumption of an "energy loss" scaling as 
E turns out to be valid for quarkonium production [43j, 
provided this "energy loss" is correctly interpreted as the 
radiated energy associated to the hard process, and not 
as the energy loss of independent incoming and outgoing 
color charges. The medium-induced radiated energy is 
parametrically similar to the Bethe-Heitler energy loss of 
a charge of mass M created in the far past, in particular 
it arises from large formation times, and the Brodsky- 
Hoyer bound does not apply in this case. 

One might wonder why large formation times tf ^ L 
do not cancel out in the medium-induced radiation, as 



in the case of a parton created in a medium (see (1.2|) 



Radiation associated with tf ^ L does not resolve the 
medium, but this does not always imply that the radi- 
ation is independent of the medium size and properties. 
For a parton created in the medium, it does. Indeed, 
in this case, radiation with tf ^ L solely arises from 
late emission off the final charged particle, and is the 
same with or without medium. This is why formation 
times tf ^ L cancel out in the medium-induced en- 
ergy loss of a parton created in the medium. In the case 
of laige-Xp quarkonium hadroproduction (and similarly, 
in the case of Bethe-Heitler radiation off an asymptotic 
charge) , although radiation with tf ^ L does not resolve 
the medium and sees it as a pointlike object, it however 
depends on the medium size via the transverse momen- 
tum broadening Aq^ cx L. Thus, the contribution with 
tf ^ L does not cancel in the medium- induced radia- 
tion associated to large-x^ quarkonium hadroproduction, 
contrary to what is assumed in Ref. [20] . 

The effect discussed here differs from that studied in 
Refs. [35]. In those studies, the nuclear suppression at 
large is explained from simple energy conservation 
arguments. In the limit — > 1, the production of ad- 
ditional particles is forbidden, and each rescattering of 
the (initial and/or final) energetic parton is associated 
to a Sudakov factor ~ {1 — Xp). Since there are more 
rescatterings in a nucleus than in a proton target, this 
naturally leads to nuclear suppression at large x^. This 
effect appears for all processes, including those (like DIS 
and DY production) where no energetic color charge is 
present in the initial or final state of the hard subprocess. 
It is argued in Refs. [22] that the effect acts similarly to 
an effective parton energy loss scaling as the parton en- 
ergy, independently of the process under consideration. 
This is clearly different from the actual parton energy loss 
(and more generally associated radiation) studied in the 
present paper. Although energy conservation is obviously 
more and more important when Xp approaches unity and 
has to be implemented in any realistic phenomenological 
model, we expect our effect to play a crucial role even 
far away from the edge of phase-space and therefore in- 
dependently of the constraints from energy conservation. 
The implications of our results on phenomenology will be 
studied in a future work. 

The paper is organized as follows. In section |ll] we 
discuss the radiated energy associated to a hard process, 
in a simple QED model. In section [II B| we consider the 
case of large angle scattering, where the radiated energy 
can be identified with the energy loss of a given charge, 
as commonly assumed. In section [il C[ we present a sim- 
ple, small-angle scattering process where this identifica- 
tion is not possible due to the presence of interference 
between initial and final state radiation. This situation 
is generalized in section |III| to the QCD case of large- 
Xp quarkonium hadroproduction, where we obtain the 
non-abelian analog of (1.5) for the medium-induced ra- 
diation spectrum. Some more detailed questions, such 
as the dependence of our results on the precise quarko- 
nium production mechanism, as well as the comparison 
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with other processes (e.g. open heavy flavour and light 



hadron production) are addressed in section IV 



To simplify the discussion in sections II B and II C we 



consider QED processes. The general conclusions drawn 
also hold in QCD. 



II. RADIATION SPECTRUM ASSOCIATED TO 
A HARD PROCESS 



B. Large angle scattering 



A. Preliminary considerations 

In the following discussion we shall consider the 
medium-induced gluon radiation spectrum, ujdi/duj, aris- 
ing from the transverse momentum broadening of par- 
tons propagating through a QCD medium. The natu- 
ral Lorentz frame to study this problem is therefore the 
frame in which the medium is static, denoted by (S). This 
frame is different whether parton propagation occurs in 
hot quark-gluon plasma (A-A collisions) or cold nuclear 
matter (say, p-A collisions). In the former case, the frame 
(S) is the center-of-mass frame of the heavy-ion reaction, 
neglecting for simplicity the longitudinal and transverse 
expansion of the plasma. The frame (S) of the latter is 
the target nucleus rest frame, boosted by a Lorentz factor 
7 — y/sj (2mp) » 1 with respect to the cm. frame of the 
p-A collision. In the static frame (S), two kinematical 
situations can be considered: 



(i) largt angle scattering, discussed in II B It corre- 



sponds typically to the production of large-p^^ par- 
ticles produced around mid-rapidity and propagat- 
ing in QGP, as illustrated in Fig. [2^; 



(ii) small angle scattering, discussed in II C which oc- 
curs in two distinct cases, namely parton propa- 
gation in cold nuclear matter (since all momenta 
are mostly longitudinal in the target nucleus rest 
frame from the large boost), see Fig.[2]D, as well as 
large-rapidity particle propagation in QGP, Fig. [2];. 




(a) 




(b) 



(c) 



FIG. 2: Large angle: (a) Particle produced at mid-rapidity 
in A-A collisions. Small angle: (b) particle produced in p~ 
A collisions, (c) particle produced at large rapidity in A-A 
collisions. 



Here we discuss large angle scattering in the frame (S), 
and first consider "p-p" collisions. The hard partonic 
process is modelled by an incoming electron, of mass m 
and energy E ^ m, scattering at large angle (90° in the 
frame (S)). The photon radiation spectrum associated to 
the hard process is obtained by calculating the photon 
emission amplitudes off the incoming and outgoing elec- 
tron lines represented in Fig [Sj 



(P) 



e~(p') 



(P) 



e" (p') 



FIG. 3: Model for large angle scattering in QED in "p-p" 
collisions. The blob represents the hard process, the photon is 
radiated by the incoming (left) and outgoing (right) electron. 

In the limit of soft photon energy uj <^ E, the photon 
emission vertices factorize from the hard amplitude, and 
one can write the photon radiation intensity as 

2 



i=l,2 



• k 



P_ 
P' 



(27r)32a; 



(2.1) 



where p — {E,p) and p' — [E' ,p') are the four-momenta 
of the incoming and outgoing electrons, k — (w, fc) is 
the photon four-momentum, Si = (0,ei) are two trans- 
verse photon polarization vectors, and e is the QED cou- 
pling. The angular integral of (2.1 1 can be performed 



exactly [23], yielding the soft photon energy spectrum 
uj^ = — [n{v) + n(v') +I(v,v' ,v ■ v')] , (2.2) 

where 7?.(f) (respectively TZ{v')) corresponds to the 
square of the emission amplitude off the initial (final) 
electron line, and I{v,v',v ■ v') stands for the interfer- 
ence term. It is worth recalling that these three terms 
only depend on the initial and final electron velocities 
V = p/E and v' = p'/E'. We have 

1 

2v 



n{v) = 



In 



T{v,v' 



l + v 
1 - V 

n{v) 



- 1 



(2.3) 
(2.4) 



where /3 is the relative velocity of the final electron in the 
rest frame of the initial one. 



/?(«, v' , V ■ v') 



(l-v2)(l-u'2) 



{1-v-v') 



l\2 



(2.5) 
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Hence, the energy spectrum reads 



dl _ 2a 
du) n 



7^(/3) 



(2.6) 



This has the following simple interpretation. In the rest 
frame of the initial electron, the final electron is suddenly 
accelerated (to the velocity and the spectrum must 
be fully given by the square of the emission amplitude off 
the final electron line. 

We focus on the ultrarelativistic limit where v, v' ^ 1. 
In this limit, 



n(v) 



1 



In 



1 



1 



oo . 



(2.7) 



showing that the squares of the initial and final emission 
amplitudes suffer from a logarithmic coUinear singular- 
ity. This corresponds to DGLAP radiation ^2T, which in 
QCD is resummed in the initial (final) parton distribu- 
tion (fragmentation) function. When v, u' — >■ 1 and the 
angle between v and v' is fixed, we have /? — t- 1, and the 
interference term reads 



X(z;, w', -I 



1 + ln 



1-v 



(2., 



This shows no logarithmic enhancement when v and v' 
are quasi-orthogonal (u • i/' <C 1). 

Thus, in the kinematical situation of Fig. [Sj the asso- 
ciated radiation spectrum is dominated by the squares 
of the initial and final emission vertices, due to a loga- 
rithmic enhancement of coUinear emission. At ultrarel- 
ativistic energies, the radiation spectrum associated to 
large angle scattering is given, in the leading logarith- 
mic approximation, by the sum of emission spectra of 
well-identified "partons" , namely those of an incoming 
(suddenly decelerated) and outgoing (suddenly acceler- 
ated) electron. We stress that the dominance of coUinear 
logarithms implies the dominance of large photon forma- 
tion times i/ in the radiation spectrum. For instance, 
the logarithm in (2.7) arises from the angular domain 
0^-^ <^ 9^ <^ 1, where 9 is the photon emission angle with 
respect to the initial electron, and 6^ = rn^ /E'^ = 1 — u^. 
The photon formation time tf ^ l/(wf?^) is thus very 
large in the ultrarelativistic limit. 

The identification between radiation associated to 
large angle scattering and parton radiation trivially ex- 
tends to the case of heavy-ion collisions. Consider some 
(abelian) QGP of finite size L produced in a central A-A 
collision, shortly after the time of the hard "partonic" 
process of Fig. [3j In the hot environment, the final elec- 
tron undergoes soft rescatterings, which will modify the 
radiation spectrum associated to the hard process. Let us 
assume that large photon formation times tf ^ L domi- 
nate in the radiation spectrum. Then, photon radiation 
does not probe the medium, and the radiation amplitude 
is given by the sum of the emissions off the initial and 
final electron lines, as in the absence of a medium. The 
in-medium rescatterings affect dominantly the direction 



of the final electron velocity, and the radiation spectrum 
in the presence of a medium is given by 



dl 

duj 



mcd 



2a 

TT 

2a 

TT 



[n{v) + n{v') + i{v, v',v-v' + Sv- v')] 



[n{v) + n{v')] , 



(2.9) 



where dv-v' arises from the in-medium modification of 
the final electron direction. For v-v' = and dv-v' <C 1, 
the interference term has no logarithmic enhancement, 
and arises from small formation times. Hence, it must 



be dropped in the first line of (2.9), which was obtained 



assuming tf ^ L. As in the vacuum case, the spectrum 
is dominated by the coUinear logarithms of the squared 
terms, and thus by large photon formation times tf ^ L. 

In the phenomenological analyses of nuclear effects, 
production rates in A-A collisions are normalized by the 
similar rates in p-p collisions. Then, jet-quenching does 
not depend on the total radiation associated to the hard 
process, but rather on the additional radiation occurring 
with a medium, when compared to the "vacuum" , p- 
p case. In the following we will focus on the so-called 
medium-induced radiation spectrum, obtained by sub- 
tracting the vacuum contribution from the in-medium 
spectrum. Subtracting (2.2) from (2.9), the dominant 



(coUinear) terms cancel out, removing the contribution 
of large formation times tf ^ L in the induced spectrum. 
Hence, the medium-induced radiation associated to large 
angle scattering must originate from photons with limi- 
ted formation time, tf < L. Since the expression (2.9) 



was obtained assuming that large formation times domi- 
nate, it is not adequate to derive the medium-induced 
spectrum. More work is needed to correctly derive the 
latter, as we now recall. 

Photons with tf < L can probe the medium size L, 
and the radiation amplitude off the scattered electron is 
not simply given by the coupling to the external lines. In 
general the induced spectrum depends on the details of 
the electron rescatterings in the medium. For the pur- 
pose of the present discussion it is sufficient to consider 
the case of a small medium of size L ^ A, with A the 
electron mean free path in the medium, so that the scat- 
tered electron undergoes at most one elastic scattering. 
The radiation amplitude induced by such a scattering is 
given by three diagrams represented in Fig. |4j 



FIG. 4: Photon emission diagrams in "A-A" collisions, where 
the outgoing electron rescatters in the medium. 

The diagram where the photon is emitted from the 
internal electron line {i.e., between the hard scattering 
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and the soft in-medium rescattering) , is negligible when 
tf ^ L and thus in the total spectrum (2.9), but becomes 
important when tf < L. Squaring the sum of the ampli- 
tude of Fig. [4j and subtracting the vacuum contribution 
(see Fig. [3|, we find that the square of the emission am- 
plitude off the incoming electron cancels out, and that 
the (induced) interference between the emissions before 
and after the hard production vertex is negligible. Thus, 
in order to derive the medium-induced radiation spec- 
trum associated to large angle scattering, we can discard 
the incoming electron, and simply consider the radiation 
amplitudes of Fig. |5] (for p-p collisions) and of Fig. [6] 
(for A-A collisions), where the blob represents the hard 
scattering. The medium-induced spectrum is that of an 
electron "created" at the hard scattering time t = 0, and 
propagating through the medium. 



production 
point 



term c>c 9 corresponds to the graph in Fig.|6^ and the term 
(X 9' to the graph in Fig. ^p. 9' is the angle between the 
photon and the electron direction after the rescattering, 
9' = 9 — 9s, where 9s = (1±/E is the electron scatter- 
ing angle {9s ^ 1) and q± the transverse momentum 
exchange in the elastic scattering. 

The medium-induced spectrum is obtained by sub- 
tracting the vacuum from the in-medium spectrum, and 
multiplying by the single scattering probability ^ L/\, 



duj 



L 
A 



a 



72 _ 72 

^ mod vac 



(2.12) 



Using the expressions of Jvac and Jmcd, and the identity 









2 




2' 


1 d'9 




9' 

-Qi2 




9 

02 





-0, 



the induced spectrum can be expressed as 



(2.13) 



FIG. 5: An electron produced in the vacuum at an initial time 
t = by a hard process radiates a photon. 



(a) 



(b) 



FIG. 6: Photon radiated by an electron produced at t = in 
the medium. 

In order to display the parametric dependence of the 
medium-induced spectrum, let us review the derivation 
of Ref. (for L <^ X). In vacuum, the radiation spec- 
trum is obtained by squaring the diagram of Fig. [5] Ne- 
glecting the electron mass and working in the small angle 
approximation, the spectrum reads 



UJ 



dl 

duj 



a 



d^9 J,,^^ 



Jv, 



02 



(2.10) 



where the vector 9 = k±/uj denotes the photon "angle" 
with respect to the final electron. The radiation spec- 
trum induce d by a single in-medium rescattering is ob- 



tained from (2.10) by replacing J, 



J n- 



where the 



in-medium emission current reads [1] 



dl 

dijj 



ind 



L 2a 



1 — cos 



wLo02- 
2 

(2.14) 

The identity (2.13) states that the radiation (inte- 
grated over angles) occurring after the electron scattering 
{i.e., the square of Fig. [6]d) is identical to the vacuum 
DGLAP radiation. (This is the small angle expression 
of the cancellation of the 'R-{v') term between (2.9) and 



(2.2).) This is because the radiation after the soft rescat- 
tering depends only on the photon angle 9' = 9 — 9s 
with respect to the final electron. In the angular inte- 
gration, the shift 9 9 + 9s removes all dependence on 
the scattering, and thus on the size or properties of the 
medium. Hence, there is an exact cancellation when re- 
moving the vacuum part to obtain the medium-induced 
spectrum, corresponding to the cancellation of large for- 
mation times tf ^ L, and the induced spectrum is thus 
dominated by < L. This is a crucial point, which is 
at the basis of the drastic difference between the present 
situation and the small angle scattering case studied in 
section III CI 

To explicitly verify in the present situation that the 
medium-induced energy s pectru m is indeed dominated 
by ^ L, let us simplify (2.14) by averaging over 9s as 



in Ref. [4|. First, we average over azimuthal directions of 
9s using 



2tt \ 02 



9(0,^ - 0^) . (2.15) 



mod 



02 



1 



-iuLoS /2 



(2.11) 



Here Lq is the distance travelled by the electron between 
its production point and the scattering. When the elec- 
tron is produced in the medium, we have Lq ^ L. The 



We then average over 02 using the (screened) Coulomb 
scattering probability distribution 



(2.16) 
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where /x is the typical value of the transverse momentum 
exchange q± in Coulomb scattering. We obtain 



dl 

duj 



ind 



L 
A 



2a 



de- 



1 - cos(wLo6'V2) 



At high energy we ha ve ^"^ /E"^ ^ l/(a;Lo 



gular integral in (2.17) is saturated by 6 
by formation times 

1 



(2.17) 

), and the an- 
1/{ujLq), i.e., 



(2.18) 



The medium-induced energy spectrum reads 



dl 

dio 



ind 



2, ,2 



and integrating this spectrum up to w 
medium-induced electron energy loss 



AEu 



A 



(2.19) 

E yields the 

(2.20) 



We see that the medium-induced energy loss is energy in- 
dependent and much smaller than the "DGLAP energy 
loss" (obtained by integrating (2.101) which is propor- 
tional to E. This is a direct consequence of the constraint 
tf < L for medium-induced radiation. We stress that the 
medium-induced spectrum is well-defined in the m 
limit, and thus collinear safe. This is because the con- 
tributions to the vacuum and in-medium spectra arising 
from small angles 0^ ~ — 0, i.e., from large formation 
times, solely arise from the DGLAP terms which cancel 
in the induced spectrum. 

The main lesson of this section can be summarized 
as follows: The medium-induced spectrum associated to 
large angle scattering can be identified to that of a charged 
particle created at t — 0. Large formation times cancel 
out in the medium-induced spectrum, leaving only tf < L, 
and resulting in an energy-independent and collinear safe 
medium-induced energy loss. 



We want to study the radiation spectrum associated to 
a generic small angle scattering process, and see how the 
spectrum is modified when going from "p-p" to "p-A" 
collisions, within a simple QED model. Our arguments 
are quite general and also apply to non-abelian radiation 
(despite slight parametric differences), as we will explic- 
itly see in section |III| in the QCD case of quarkonium 
production. Our QED model for small angle scattering 
is depicted in Fig. [T] 




FIG. 7: QED model for the production of a massive particle 
(muon) at small angle. The outgoing electron and anti-muon 
are soft and do not participate to the hard scattering dynam- 



An electron of mass m and energy E ^ m undergoes a 
hard scattering, from which a muon of mass M ^ m and 
energy i^' ~ ^ Af emerges. The particles produced 
in the final state in conjunction with the energetic muon 
(which are actually required here from lepton number 
conservation) are soft. Radiation off those soft particles 
(for instance, the outgoing /x+ and e~ in Fig. [7) can be 
disregarded, and only the incoming and outgoing ener- 
getic charges are relevant for our purpose. 

The radiation spectrum associated to the process of 
Fig. [T] is easily obtained by assuming that it arises from 
photon formation times satisfying i/ 3> ^hard: which can 
be checked a posteriori The radiation amplitude is 
then given by the two diagrams of Fig. [Hj 



C. Small angle scattering 

We now discuss small angle scattering processes in the 
frame (S). In the following, the invariant mass AIx of the 
partonic final state sets the hardness of the process. We 
work in the target rest frame, where small angle scatter- 
ing is modelled as an incoming energetic charge of mo- 
mentum p = (E,0±,pz) scattering off the target with 
limited momentum exchange q. Creating the invariant 
mass Mx requires some non-zero longitudinal momen- 
tum transfer g^, 



Ml 
E 



(2.21) 



which using the uncertainty principle is related to the 
coherence time ihard of the hard partonic process. 



^hard 



1 



E 

Ml 



(2.22) 




(a) 



(b) 



FIG. 8: Photon emission by the incoming electron and out- 
going muon, in the QED model for small angle scattering. 



The radiation spectrum is thus given by the expres- 
sion (2.2) (or (2.6)), where v' is now the velocity of the 
outgoing muon. For small angle scattering, v ■ v' ~ 1. 
Denoting Og -^1 the angle between v and w', and using 
= in^/E^ = l-v^ < 1 and_^ = M"^ / E^ ~ 1 



1, the relative velocity /3 (see (2.5)) between the electron 
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and muon can be approximated by 




(2.23) 



Since M ^ m, we have (3^1. Using (2.8) we find 



ti', V ■ v') ~ In (1 



In 



^^1 



Hi 



(2.24) 

Thus, the interference of the two diagrams of Fig. |8] is 
enhanced by a colhnear logarithm, in contrast with the 
large angle scattering case discussed in sect ion |II B| 
The radiation spectrum is obtained from (2.2 1 by using 



(2.7) and (2.24). For p-p collisions. 



dl 

doj 




, (2.25) 



The first two terms stand for the DGLAP radiation of the 
electron and muon respectively, and the third term is the 
interference. Note that the interference term largely com- 
pensates t he D GLAP terms. We easily verify that the 
spectrum ( 2.25 ) arises from formation times tf ^ E/M]_ , 
i.e., t f ^ ihardj as initially assumed. We stress that 
(2.25) is the spectrum associated to the hard process. 



Here, contrary to section II B[ it cannot be interpreted as 
the radiation spectrum of a well-identified charge, due to 
the presence of a large interference term. 

Inp-A collisions, the angle between the outgoing muon 
and incoming electron tends to increase, due to trans- 
verse momentum broadening in the target nucleus. For a 
nuclear target of sufficiently large size L ^ X (with A the 
electron mean free path in the target), transverse momen- 
tum broadening is given by the random walk estimate 
~ where /i is the typical transverse ex- 

change in a single scattering. Thus, ^'slpA — ^s|pp + ^^sj 
where A6'^ = Aq'^/E'^ is the angular broadening. For- 
mation times tf ^ thard 3> L being dominant (as in 
the p-p scattering case), the spectrum is equivalent to 
that associated to a single effective scattering of trans- 
verse exchange = g^.lpp + ^'zl- We focus on the limit 



A(72^ ^ ^ilpp' ^^^^ '^^^ hard production process is 
not affected by soft rescatterings. The broadening of Bg 
induces a slight modification of the interference term in 
p-A collisions, as compared to p-p (see ( 2.24[ )), leading 
to the spectrum 



dl 

dio 



a 








m 











-2 In 



92 



-i IpA 



, (2.26) 



with again a large contribution from the interference 
term. 



The medium-induced spectrum is obtained by subtract- 



ing (2.25) from (2.26) 



dw 



ind 

2a 

TT 



2a 



In 1 



In 1 



Ml 



92 



pp y 
2 



2a Ag2 



(2.27) 
(2.28) 



The expression (2.28) was obtained using (2.2), which 
arises from an exact angular integration. At high en- 
ergy, all relevant angles {0,9s, 9m, Om) are small, and one 
could have as well worked in the small angle approxima- 
tion from the beginning. In view of our generalization 



to the non-abelian situation in section III it is useful to 



mention how the expression (2.27) arises within this ap- 
proximation. Expressing the square appearing in (2.1 ) in 
the small angle limit, it is easy to check that the radiation 
spectrum as sociat ed to our small angle scattering process 
is given by (2.101, with Jvac replaced by the in-medium 



"emission current" 



J n- 



QI2 



92 



(2.29) 



The second term arises from Fig. [§1 ( emission off the 
electron line) and the first, where 0p^ — 9 — Og with 
9^s — IpA' fro™ Fig. jsjj (emission off the muon line). 
The medium-induced spectrum thus reads 



d]_ 
duj 



ind 



a 



72 _ j2 



2a 



d^9 



92 



Q/2 



(vac) 



(2.30) 
(2.31) 



where we set m 
obta: 
have 



obtained by replacing 9^ 



and the vacuum contribution is 

A ^ ^'1 ■ 

pA ^ I pp 



From (2.15) we 



27r 6»'2 



e(6'2 - 9 



^) 

s J J 



and we arrive at 
dl 

did 



ind 



2a 

TT ./6)2| 



IpA d9 

W+9i 



(2.32) 



(2.33) 



M 



from which the expression (2.271 follows 



The medium-induced spectrum (2.33) is dominated by 
formation times 



t 



1 



1 



/ 



(2.34) 

Hence, 3> thard 3> L, justifying our initial assumption. 
Contrary to the case studied in section [II B[ large photon 
formation times do not cancel out and dominate in the 
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medium-induced spectrum. This justifies approximating 
the radiation spectrum as that associated to a single ef- 
fective scattering, and using (2.2 1 throughout the discus- 
sion. The induced spectrum arises solely from the (in- 
duced) interference term, and cannot be identified with 
the induced spectrum of a given charged particle, either 
the electron or the muon. To stress this point, let us 
note that if the incoming (outgoing) energetic particle 
carries a charge e (e') in units of the electron charge, the 
medium-induced spectrum will be given by (2.27) multi- 
plied by ee'. 



to be contrasted with the result (2.20) for large angle 



The induced spectrum (2.27) is very similar to the 



Bethe-Heitler (BH) radiation spectrum off an asymptotic 
muon crossing a nucleus. The latter is obtained from 
(pjl) by setting m = M in (12. 23^. For 6*^ = 0, we have 



scattering. We note that the parametric dependence of 
the radiated energy (2.38) is formally the same as that of 



the energy loss of an asymptotic charge of mass M± , see 
( p36l ). 

To summarize this section: In the case of small-angle 
scattering, the medium-induced photon radiation spec- 
trum associated to the hard process arises from the inter- 
ference between initial and final state radiation. As such, 
it cannot be identified with the medium-induced spec- 
trum of well-defined charged particles. Nevertheless, the 
spectrum is qualitatively the same as that of an asymp- 
totic massive charge crossing the medium. In particular, 
it arises from formation times tf ^ L, resulting in a 
medium-induced radiated energy scaling as E . 



/3 = and the BH spectrum vanishes. For Qg ^ 0, 



1 



(l + a/2)2 



92 



(2.35) 



and the function TZ{fi) appearing in (2.6) and defined 
by (2.3) is very well approximated by In (1 -t- a/3), which 



has the same limiting behaviours as 7?,(/3) when a ^ 1 
and a 3> 1. The soft photon radiation spectrum off an 
asymptotic massive charge thus reads [1] 

3^ 



dl 



duj 



2a 



In 1 



sell 



(2.36) 



In the small angle approximation, the result (2.36) would 



be obtained from (2.29) and (2.30) by setting M — m, 

2 




(2.37) 



Performing the angular integral one recovers (2.36) 



The BH spectrum ( 2.36 ) is qualitatively similar to the 
radiation spectrum (2.27) associated to the hard produc- 
tion process of Fig. 7] In particular it is dominated by 
large formation times tf ^ L. It is worth noting that 
the radiation spectrum of an asymptotic charge can be 
obtained from the induced spectrum (2.12), (2.14) of a 



charge produced in a medium by moving the production 
point to the far past, Lq ^ 00. This is another way to 
understand why large formation times contribute in the 
present case. The typical formation time contributing to 
the induced loss of a suddenly accelerated particle (see 

When Lq increases at fixed 



IIB) is tf ^ La^ L. 



section 

L,tf increases to larger and larger values, and eventually 



saturates a.ttf'^ 1I{ijj6\j). 



Integrating the energy spectrum (2.28) up to a; ^ i? 



we obtain the medium-induced radiated energy (rather 
than energy loss) associated to the small angle scattering 
process. 



ind 



Ml 



I4?_ 
XMl 



E. 



(2.38) 



III. LARGE QUARKONIUM PRODUCTION 

In this section we consider the QCD process of quarko- 
nium (hadro)production at large- a;^, in p-p and p-A col- 
lisions. We discuss the medium- induced radiation spec- 
trum associated to the hard process, given by the differ- 
ence between the soft gluon radiation spectra associated 
to vacuum {p~p) and in-mcdium (p-A) production. We 
will show that the features of gluon radiation are simi- 



lar to those obtained in the QED model of section II C 
(despite some difference in the parametric behaviour of 
the photon and gluon energy spectra). In particular, the 
medium-induced spectrum arises from large gluon forma- 
tion times, leading to a medium-induced radiated energy 
scaling as the quarkonium energy. 



A. Model for large-a;^ quarkonium 
hadroproduction 

In order to single out the main features of large- a;^ 
quarkonium production, we use several simplifying as- 
sumptions. As mentioned in the Introduction, we focus 
on the partonic process gg — >■ QQ, where the final quark 
and antiquark momenta are similar and quasi-collinear. 
At large xp = xi — X2 — x\ and in the target rest 
frame, the incoming gluon splits into the QQ pair, to 
which it transfers most of its energy E, equally shared 
between the quark and antiquark. The coherence time 
of the hard process is given by (2.22) with Mx = 



Afj^ ^ Aqcd denoting the quarkonium transverse mass. 
Thus, ihard ~ E / M\, corresponding to the g — >• QQ fiuc- 
tuation time in the target rest frame. 

In the following we will not have to specify the pre- 
cise mechanism for the transition between the QQ pair 
and the quarkonium bound state. We will however as- 
sume that the QQ pair remains a color octet for a time 
ioctet ^ E/M^. This is a reasonable assumption for 
the Color Octet Mechanism (COM) 25j, or in the Color 
Evaporation Model (CEM) [26 for quarkonium produc- 
tion. In the Color Singlet Model (CSM) [23, the quan- 
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turn numbers of the quarkonium bound state (in par- 
ticular its color neutrality) are fixed in the perturba- 



tive process, i.e., 



CSM 
octet 



t 



hard 



E/M[ 



'j_ instead of 



^octot ^ E/M'j_. Our results might nevertheless also have 
implications on quarkonium production within the CSM, 
see section [lVC|for a discussion of this point. 

In section IIIB| below, we will study the medium- 
induced gluon radiation spectrum associated to large- a;^ 
quarkonium hadroproduction in the u domain where 



E 

Ml 



< t 



f 



< to 



(3.1) 



the first inequality implying that radiation does not 
probe the hard gg — > QQ partonic process, and the se- 
cond that the QQ pair remains color octet during the 
overall radiation process. We will also assume that the 
hard scale Mj_ = ^ + is the largest scale after 
the quarkonium energy E, where q± is the transverse 
momentum of the QQ pair. In particular, the radiated 
gluon has k± <C M±,q±. Under those conditions, the 
radiation does not affect the kinematics of the hard pro- 
cess, and does not probe the transverse size ^ O (1/M_l) 
of the QQ pair |45j, and the latter thus appears as a 
compact colov octet state, i.e., as a (massive) gluon [55]. 

For our purpose, the quarkonium production process 
is effectively equivalent to the process depicted in Fig. |9j 
where an energetic "massive gluon" is produced in the 
hard scattering of an incoming (massless) gluon. We 
stress that the incoming gluon should not be coherent 
with the other partons of the projectile proton for our 
perturbative treatment to be meaningful. This requires 
the typical transverse momenta to probe the proton size, 
fc_L,g_L > Aqcd- 



iiiiiiiiiiiiiiMililllilllillll 



X 



arises from large (gluon) formation times, see (3.11, in 
particular tf ^ L at large E. Under those conditions, 
gluon radiation does not probe the hard process, and only 
emission vertices before or after the hard blob of Fig. [9] 
have to be considered. Similarly to Ref. |3], we use semi- 
heuristic arguments to derive the radiation spectrum. 

The emission vertex off the final gluon line, obtained by 
factoring out the amplitude without radiation (see Fig. [9| 
from the amplitude with radiation, has the same struc- 
ture as in QED (up to some implicit color factor). 



nal 



0,2 
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<1± 

E 



(3.2) 



The term 0\j in the denominator arises from the fact that 
the pointlike color octet QQ pair (i.e., the final "gluon" 
in Fig. |9| has a mass M, and 9' is the emission angle 
with respect to the final "gluon" direction {9s — q± /E is 
the energetic gluon scattering angle). 

The emission vertex off the initial (massless) gluon line 
is a priori more difficult to extract, since the radiated 
gluon can rescatter in the target. Let us consider the 
process in a p-A collision. A generic diagram for the 
radiation amplitude is shown in Fig. |10[ Compared to 
the process in p-p scattering, the radiated gluon suffers 
some transverse momentum broadening Ag^ ^ (L/X) /i^ 
|46| . The emission vertex can be obtained heuristically by 
noting that just after the emission, the angle between the 
radiated and initial gluon momenta is 9 " = 9 — 9g, where 
9g = A(fj^/a; is the rescattering angle of the radiated 
gluon in the target The fact that in Fig. 10 the 



radiated gluon might couple to the rescattering (Agj^) 
gluon but not to the hard {q±_) exchanged gluon arises 
from our assumption of soft gluon radiation, which should 
not affect the kinematics of the hard production process, 
and thus factorize from it. Clearly, our picture requires 
k±,Aq± <C q± to be vaUd. Thus, effectively. 



litial 



UJ 



(3.3) 



FIG. 9: QCD model for quarkonium hadroproduction at large 
. The "massive gluon" turns into a color-neutral object 
(and eventually, into a quarkonium) on longer time scales in 
the GEM and GOM. 



B. Medium-induced radiation spectrum 

We now want to derive the medium-induced radiation 
spectrum associated to large-x^^ quarkonium production, 
focusing on the case of a medium of size L 3> A, where 
A is the gluon mean free path between successive elastic 
scatterings in the target nucleus. As in the QED case 
studied in section II C[ we anticipate that the spectrum 




p'=M^ 



^QQQmmm 



OOOQOOOQO, 



XX 



FIG. 10: Generic amplitude for gluon radiation in the model 
of Fig. [9] In QCD the radiated gluon also rescatters in the 
medium. 
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The medium-induced spectrum associated to large-x^ 



quarkonium production is obtained by multiplying (3.2| 



and (3.3) [JT], integrating over cP9, and then subtracting 
the "vacuum" contribution (formally obtained by Aq± — >■ 
0_l) corresponding to p-p scattering. After a simple shift 
in the angular integration variable 9 we obtain: 



cOq / to = 0.1 
cOf, / to = 0.3 



dl 

dui 



ind 



J e^ + ol 



M 



— (vac) 



(3.4) 



The latter expression is similar to the QED expr ession 
(2.31), up to the replacement 9s — > 9s ~9g. Using (2.32|, 
the angular integration gives: 



dl 

duj 



^ )1 PA d9^ 

ind J (e,-e,Y\ ^ -^"M 

ipp 

^ ^lnfl + ^ 



2^n2- (3-5) 



(3.6) 



In order to obtain the second line, we approximated 



pA 



PP 



IpA 



(3.7) 



Ipp £;2 



Comparing to the QED result (2.28), we observe that 



the QCD spectrum (3.6) involves a new scale. 



dio 



Nras 



In 1 



ind 



(3.9) 
(3.10) 



above which the spectrum is suppressed as ~ 1/w^. 

The medium-induced radiated energy is obtained by 
integrating the spectrum (3.9 1 over w, up to ^ E. 



However, due to fast convergence for a; > a), the integral 
is well approximated by replacing E by infinity. We find: 



A£:| 



ind, large xp 



NrUs 



Ml 



E. (3.11) 



The medium-induced radiated energy arises from giuon 
energies uj ^ Co E, and thus from gluon formation 
times tf '-^ E'^/{u;M]_) E/M]_, as initially assumed. 
Note also that fc_L ~ ~ ujM_i/E <C M± and thus 
the radiated gluon cannot probe the transverse size ~ 
O (1/M_l) of the heavy QQ pair. 



The medium-induced radiated energy (3.11) scales as 



the quarkonium energy, due to a dominant contribution 
from large gluon formation times. Similarly to the QED 
case, it cannot be attributed to any well-identified par- 
ton, since it arises from the interference between emission 
vertices off an incoming massless gluon and an outgoing 




(O / to 



FIG. 11: Medium-induced soft gluon radiation spectrum 
(3.13 1 associated to large quarkonium production, for var- 
ious values of luq/oj. We chose N^cXs/t^ ~ 1/2 to draw the 
figure. 



massive pointlike color octet. However, (3.11) exhibits 



the same parametric dependence as the radiative loss of 
an asymptotic color charge of "mass" M± undergoing a 
single effective scattering of momentum transfer Ag^ [3] . 



(3.8) jjj particular, it behaves as 1/M± rather than 1/M 



in QED (see (|2^). 

Finally, in view of phenomenological applications, let 



us determine the validity range of the spectrum (3.9 1 



First, the spectrum was derived in the soft giuon ap- 
proximation, Lo <^ E. Second, and most importantly, in 
QCD our perturbative derivation is meaningful provided 
fcj^ > Aqcd, as already mentioned in section HI A Us- 
ing fcj^ ~ uj9 and the approximations (3.7), (3.8), the 



spectrum (3.5 1 is reexpressed as 



LjJ 



dl 

duj 



Nras 



ind JMaxix^q'^^.A'Qco 

where the constraint k± > A 
into account, and x = 



dkl 



kl 



(3.12) 



Lj/E. 



QCD is now taken 
Note that Aqj_ 



{L/X)fi ^ ^QCD fo'^ large nuclei. Approximating 
dl 



duj 



Nctts , /I + w^/w^ 



In 



TT 



1 + a;2/w2 



^ Aqcd 



Wo 



E <Lh . 



(3.13) 
(3.14) 



The spectrum is shown in Fig. [TT] from which we see 
that Eq. (3.9) is valid down to a; ~ ojq, below which 



it becomes smoother and the logarithm saturates at 
In {u^^Ju^D = In (Aqi/A^jcD) ^ In [L/X). 

With the constraint k±^ > Aqcd , the average medium- 
induced radiated energy now reads 



ASI; 



ind, large Xf 



NcOs (a) - cjo) ■ 



(3.15) 
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The result (3.11) thus receives a reduction factor 



(3.16) 



which is about 0.5 when L/X — A and /i ~ Aqcd- 



IV. COMPARING DIFFERENT PROCESSES 

The main message of our study can be stated as fol- 
lows. In a hard process involving incoming and outgo- 
ing energetic charges (which do not have to be identical) 
being quasi-collinear in the rest frame of the medium, 
the associated, medium-induced radiated energy - ob- 
tained by subtracting the contribution in p-p from that 
in p-A scattering - arises from the interference be- 
tween the emission amplitudes off the initial and fi- 
nal charges, and is dominated by large gluon formation 
times. The medium-induced radiated energy is similar to 
the medium-induced energy loss of an asymptotic charge, 
in particular it scales as the energy, which might have im- 
portant consequences on the phenomenology. 

In this section we discuss various hard processes off 
nuclear targets and examine for each of them whether 
the medium-induced energy loss is similar to that of 
an asymptotic charge or not. The quantitative effects 
of medium-induced gluon radiation on large-x^ quarko- 
nium, open heavy flavour, and light hadron hadropro- 
duction will be addressed in a future study. Since the 
gluon radiation spectra associated to these processes are 
similar (they only differ by the precise value of the scales 
uj and luq appearing in (3.13l), we anticipate that similar 
(and large) nuclear suppressions will be obtained in those 
cases. 



A. DIS and photoproduction 

In DIS on nuclei, the medium-induced energy loss is 
that of the parton knocked by the virtual photon. In 
this case, the expected energy loss is that of a charged 
particle produced at t = 0, as assumed in various phe- 
nomenological analyses p9] . 

The situation is a little more complex in photopro- 
duction reactions. As long as the initial photon directly 
couples to the quark in the target nucleus (at leading 
order), the process resembles DIS and the energy loss is 
that of the struck quark. In photoproduction however, 
the hadronic structure of the photon can also be resolved: 
in this case a colored parton stemming from the photon 
participates to the hard scattering dynamics, making re- 
solved photoproduction very similar to hadroproduction. 
Since the time scale to resolve the parton inside the pho- 
ton is typically much larger than the time at which the 
hard process occurs, one could expect an interference be- 
tween the emission amplitudes off the initial and final- 
state parton, leading to a medium-induced energy loss 



(or, more accurately, radiated energy) scaling like the 
parton energy. In this respect, it would be extremely 
valuable to investigate the nuclear dependence of forward 
jet and hadron production in direct vs. resolved photo- 
production in nuclei. 



B. Drell-Yan production 

As already mentioned, in the absence of color charge in 
the partonic final state, the energy loss in Drell-Yan pro- 
duction is expected to be that of a suddenly decelerated 
parton, i.e., independent of its energy |48j . as is the case 
in large angle scattering (section II B I . The consequence 
is that the effects of energy loss should play almost no 
role in Drell-Yan production in high-energy p-A colli- 
sions, unlike what is assumed in the model by Gavin and 
Milana [19 . 

The energy loss scaling as E would come into play only 
if another energetic charged particle is produced in the 
final-state, in association with the virtual photon. Such 
a situation occurs in DY-|-jet production in p-A or A-A 
collisions, with the jet produced at large rapidity. 



C. Heavy- flavour hadroproduction 

In the case of large-x^ quarkonium hadroproduction, 
our result (3.11 ) holds provided gluon radiation has time 



to be formed before the QQ pair turns color singlet, see 
As we mentioned in section 



(3.1) 



III A this assump- 



tion is justified in the GEM and GOM (provided color 
octet contributions to quarkonium production dominate 
in the latter case). As a consequence, one could expect 
in this case a different nuclear dependence of quarkonium 
hadroproduction at large Xp and quarkonium photopro- 



duction at large z (see section IV A) 



However, we should recall that the precise dynamics 
of quarkonium production is still unknown. Indeed, no 
proposed model can explain all features of the data on 
quarkonium production. Although the GSM alone seems 
to be ruled out by hadroproduction data, it is worth keep- 
ing in mind that color singlet contributions might not be 
negligible in some kinematical regions |30j . 

We have to stress that our main result (|3.9|) (or more 



accurately (3.13 1) does not apply to such color singlet 
contributions, for which toctct ~ ^hard ^ E/M'^, in con- 



trast to our assumption (3.1). In the GSM, the outgo- 



ing color octet QQ pair is too short-lived to allow for 
our interference contribution to gluon radiation. How- 
ever, depending on the quarkonium quantum numbers, 
the GSM may require the quarkonium bound state to 
be produced in conjunction with a hard gluon, as for 
instance in J/ip production. In this case, the hard pro- 
cess looks like small angle scattering of an energetic color 
charge - even though the final charge (the hard gluon) 
is distinct from the triggered QQ pair. The medium- 
induced radiation spectrum (3.13), with E interpreted 
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as the energy of the hard gluon (and M set to zero in the 
expression of M±), might thus indirectly affect the J /ip 
production rate. In other cases, hke Xc2 production in 
the CSM, where no associated hard gluon is required (at 
leading order in as), the radiation spectrum (3.13) will 
not apply. 

In the COM the Xc2 state - which radiative decays 
contribute quite significantly to J/ ij^ production ~ is pro- 
duced predominantly as a col or sin glet object. Therefore, 
we do not expect our result (3.11 ) to apply to this state, 
leading to a much smaller Xc2 suppression (as compared 
to J/V") at large x^. This prediction, which could be 
tested at fixed-target experiments at Fermilab [3T], is at 
variance with other nuclear effects such as nuclear PDF 
effects or intrinsic charm which do not depend on the 
identity of the final charmonium state [16] . Another pre- 
diction can be made regarding the T suppression. From 
the mean radiated energy (3.111, /S.E ~ 1/M±, the sup- 
pression of T production, for which no measurements 
have been performed yet at large Xp , is expected to be 
less pronounced than that of J/V"- 

We expect the medium-induced radiation spectrum 



(3.13) to be valid in open heavy flavour hadroproduction 
at large x^. In particular, the fact that in this process 
the final energetic heavy quark does not have the same 
color charge as the octet QQ pair in quarkonium produc- 
tion does not affect our conclusions. As mentioned in the 
Introduction, the similar suppression of J /if) production 
and single muon (coming from £)-decays) production in 
p-A collisions as a function of rapidity [9J supports this 
picture. 



D. Light-hadron hadroproduction 



(ii) the light-hadron suppression in p-A or A-A col- 
lisions is defined as the (normalized) ratio of the 
light-hadron production rate integrated over all 
transverse momenta (with an appropriate infrared 
cut-off). In that case, Qj_|p^ = ^^-Ipp longer 
required, leading to a potentially large energy loss 
proportional to the parton energy. The induced 
gluon spectrum is then obtained by setting M = 0, 
i.e., A/j^ <li_ in (3.13). 



In the case of the BRAHMS experiment at RHIC [TT] . 
the quenching factor for all charged hadrons in d-Au 
compared to p-p collisions at ^/s = 200 GeV/nucleon 
and at forward rapidity (0 < y < 4) is measured at fixed 
hadron transverse momentum p^ , which therefore looks 
like case (i) above. However note that the final parton 
transverse momentum is smeared because of the fragmen- 
tation process in QCD. Whenever the typical dispersion 
of the parton transverse momentum is large (as compared 
to the momentum broadening Ag^), we no longer have 
'^J-lpA ~ '^J-Ipp ^"^^ energy loss might scale with the 
parton energy E cx exp y. Another reason why compar- 
ing at fixed p^"^ might nevertheless induce an energy loss 
similar to that of asymptotic partons is the binning used 
in the experiment. As long as the typical size of the 
Pj^-bin is large as compared to the nuclear broadening, 
this case is similar to an integrated distribution over the 
transverse momenta, i.e., to case (ii). 

The NA49 experiment at SPS also measured the Xp- 
differential production cross section of light hadrons in 
p-p and p-A collisions and integrated over all transverse 
momenta '10], which corresponds to the case (ii) above. 
The strong suppression observed at large x^ therefore ap- 
pears qualitatively consistent with our expectation [49]. 



Unlike DIS, Drell-Yan and heavy-flavour production, 
in the case of light hadron production at large x^ there 
is no natural hard scale, except if the final hadron trans- 
verse momentum in a p-p collision is large enough, p^'^ ^ 
<Z_l| 3> Aqcd, in which case a perturbative description 
might hold. Regarding the expected medium-induced ra- 
diation spectrum, different cases might occur depending 
on the actual definition of the observable: 



(i) the suppression is defined as the (normalized) ra- 
tio of the light hadron rate in p-A (or A-A) colli- 
sions over p-p collisions at the same, fixed <l\\^p^ = 

q]] .In this case, the interference term in the in- 
duced radiation spectrum vanishes from the actual 
definition of the observable, see the QED expres- 



sions Eqs. (1.3) and (1.4 1, or set lS.q, = in the 



gluon radiation spectrum (3.6). There will never- 



theless be some energy loss effect due to rescat- 
tering in the medium. However, in this case the 
energy loss will be dominated by the radiation of 
gluons which resolve the medium (i.e., with a for- 
mation time of the order of the medium size L) and 
will not scale as the parton energy; 



V. OUTLOOK 

As already mentioned, the notion of parton energy loss 
is not general enough to apply to the situations with 
an important interference between initial and final state 
radiation, where the associated radiated energy should 
instead be considered. Quite interestingly, when pro- 
gressively going from large Xp (or large rapidity) to the 
central rapidity region, we expect a transition between 
a regime where the (average) medium-induced radiated 
energy scales as E and a regime where it does not. It 
would be interesting to single out some production pro- 
cess where such a transition can be put in evidence. In 
this respect, the RHIC data on quarkonium suppres- 
sion in nuclear collisions appear very promising. The 
PHENIX collaboration reported a stronger J/ -0 suppres- 
sion at forward rapidity than at mid-rapidity in d-Au |32| 
and also Au-Au collisions [33], unlike the models based 
on quarkonium dissociation in QGP. Another observa- 
tion, reported by STAR, is the smaller J/ip suppression 
when going from low to high p^ [34) . Both measure- 
ments seem consistent with such a transition between 
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the two identified regimes. More generally, we mentioned 
that some apparently unrelated nuclear suppression ef- 
fects (jet-quenching in A-A collisions, open and hidden 
heavy-flavour and light-hadron production in p-A) might 
be explained in a unified framework. 

We also pointed out in the Introduction that the scal- 
ing AE cx E predicted in the present paper for J/ip for- 
ward production would lead to the scaling of J/ip sup- 
pression with xi (or Xp ) independently of the cm. energy 
•y/fi of the p-A collision, as observed experimentally. On 
the contrary, an energy loss independent of the parton 
energy - as for instance in Drell-Yan production - would 
lead to a stronger nuclear suppression at lower beam en- 
ergy, since the typical shift of the momentum fraction in 
the proton PDF would scale like Axi cx l/i?bcam oc 1/s 
causing a breakdown of xi (or x^) scaling. Unlike the 
GM model we do not expect any scaling of DY suppres- 
sion in either xi or x^ variables. Lacking precise DY 
measurements at various beam energies, this prediction 
has not been tested yet. The future data to be collected 
by the E906 experiment (E^beam = 120 GeV) 05\ at Fer- 
milab and possibly at J-PARC (£;bcam = 50 GeV) [35] 
will allow for crucial tests of the energy loss effects dis- 
cussed in this paper, both in DY and J/ip production. 

On a more theoretical ground, it would be valuable to 
study how our results are related to QCD factorization 
issues. The physics described in our study is clearly of 
coUinear nature, as can be seen for instance from the di- 
vergence of our result (3.6) (or its abelian version ( |2.28 )) 
in the formal limit Mj_ — !■ 0. However, as we already 
mentioned, our analysis requires AIj_ ^ k±, Agj^, so 



that the kinematics of the hard process is not affected by 
(soft) radiation and rescatterings. As a consequence, the 
avera ge ra diated ener gy is suppressed by a power of M± 
(see ( |2.38 1 and (3.11)), i.e., it is of higher- twist. Thus, 
our results should not conflict with leading-twist collinear 
factorization [37]. However, we do not see why our effect 
should disappear when the kinematics is extended to the 
domain where M±, A.q± and k± are of the same order. 
As already mentioned, in QED one may consider the in- 
coming and outgoing energetic particles participating to 
the hard process to carry different charges e and e'. The 
fact that the medium-induced radiation spectrum (2.28) 



is then proportional to ee' suggests that the effect can- 
not be attributed to any purely initial (parton density) 
or final (fragmentation) effect. For the same reason, that 
the effect could be part of other (more general) universal 
objects, such as parton correlation functions (i.e., fully 
unintegrated parton densities) [38], seems also unlikely. 
It has been recently argued that factorization (even in 
a generalized, AiT-dependent formulation) is most proba- 
bly violated in some cases, for instance in the production 
of laige-pT back-to-back hadrons in hadron-hadron colli- 
sions [39) . Similarly, it is plausible that hard-scattering 
factorization is truly violated in the small angle scatter- 
ing process we have considered in the present study. 
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